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J> Abstract 

(N 

^vq The Galactic secondary cosmic ray anti-proton (p) flux calculated with different propagation models is fairly consistent with data, 
and the associated propagation uncertainty is small. This is not the case for any p exotic component of the dark matter halo (see 

^\ also the companion paper; Maurin et al. 2006). Detailed propagation models are mandatory if the ultimate goal is to explain an 
excess. However, simpler and faster approximate formulae for p axe an attractive alternative to quickly check that a given dark 
matter model is not inconsistent with the p observed flux. This paper provides such formulae. In addition, they could be used 
to put constraints on new physics in this channel, where an extensive scan of a large parameter space could otherwise be quite 
(-H expensive in computer ressources. 

Oh- 
I ■ 

Q Key words: Cosmic Rays, Diffusion equation, Anti-nucici, Dark matter, Indirect detection 
^ 'PACS: 98.38.Cp, 98.35.Pr, 96.40.-z, 98.70.Sa, 96.50S-, 96.50sb, 95.30.Cq, 12.60.Jv, 95.35.-l-d 



.5^ 1. Motivation 



The first po sitive detection of anti-protons at the end of 



the seventies ( Golden et al. . 19791 : BufRngton et al. . 1981 



Bogomolov et al. . 19821 ) boosted the interest for indirect 



searches of new physics in the low energy p spectrum. After 
25 years of efforts and improvements both in measurements 
and theoretical calculations, the data at low energy (up to a 
few GeV) arc now well accounted for ( Donato et al.1 . 12001 ). 
However, at highe r energy, calculation s tend to predict less 
p than observed (jBoezio et al.l . 120011 ) but more data are 
desirable to confirm this possible trend. 

The uncertainty of the p standard secondary flux 
( Donato et al.l . [200 ih mainly comes from nuclear physics, 
as that of astrophysical origin is small — less than a few 
percents — and is expected to decrease even more as new 
B/C measurements are available. In contrast, the p ex- 
otic flux — from exotic sources following the dark matter 
distribution — is plagued by the degeneracy in the prop- 
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agation parameters, the corresponding uncertainty be- 
ing as large as two orders of magnitude at low energy 



()Donato et al.l . 120041 ). Until propagation is better under- 



stood and constrained, it may not be worth using refined — 
and sometimes time consuming — models whose result will 
anyway be crippled by these large uncertainties. This is es- 
pecially true when one has to scan as efficiently as possible 
the huge parameter space existing, for example, in SUSY 
theories. 

The goal of this paper is to provide approximate for- 
mulae for the propagated p exotic fluxes. This is done 
in the framework of the diffusion /c onvection model wit h 



constant wind alrea d y discussed iii iMaurin et al.l (120011): 
Donato et aP (|200lf l: iBarrau et all (|2002l ): IPonato et al 
(|2004[ ). These formulae are reasonably accurate, fast to 
compute and easy to implement. They have to be thought 
as an easy-to-use tool, for phenomenologists interested 
in beyond-the-standard-model theories and wishing to 
quickly check that any new physics model on the market 
does not violate the p constraint. Once the interesting 
regions of the parameter space (for the dark matter can- 
didate) are identified, more elab orate treatments should 



be implemen ted (e.g. Galprop — [Moskalenko et al 



2005 



DarkSusy— iGondolo et al.ll2004l: or lDonato et al.ll2004 ) 



The paper is organized as follows: 

- In Sec. 2, we remind the salient ingredients of the con- 
stant wind/diffusion model for secondary and primary 
exotic anti-protons, using a Bessel expansion formalism. 

- In Sec. 3, two alternative formulations of the p primary 
exotic flux are presented. 

- In Sec. 4, all the formulations are compared and their 
relative merits discussed. 

The formulae given in Sec. 3 are implemented in numerical 
routines publicly availabl(P1. 

2. 2D model with constant wind Vc 

Cosmic ray fluxes ar e determined by the tran sport equa- 



out the paper, we use the so-called thin-disk approxima- 
tion where the gas is contained in a layer of thickness 2h = 
200 pc. Cylindrical symmetry is assumed, and the diffusion 
coefficient is constant in the whole Galaxy, 

K{E) = iSKoTZ^ (7^ = pc/\Z\e is the rigidity). 

A galactic wind Vc of constant magnitude, directed out- 
wards along z, is also included. 

2.1. Diffusion equation 

Denoting the differential density as N leads to 
d 



- 2h6{z) 
--S{r,z,E) 



(1) 



The quantity Ftot = SisM "-iSM.t^.cfgi^^ is the destruction 
rate of p in the thin gaseous disk (nisM =H, He). The 
terms b{E) and c{E) correspond respectively to a drift 
term (coulomb, ionization, adiabatic losses and reacceler- 
ation) and a diffusion term (reacceleration depending on 



the Alfveni c velocity) in energy space (e.g.. iMaurin et al 



2OOIU2OO2I ). The quantity S{r, z, E) stands for the source 



term and has three contributions: secondary (spallation- 
induced), tertiary (non-anni hilating rescattering) and pri- 
mary (exotic dark matter)— IPonato et all (|200ll l2004l ). In 
the remaining of the paper, we term energy redistribution 
the effects of drift, diffusion and tertiary contribution. 

2.2. Propagation parameters 

The parameters of the model relevant to the present 
study are i) the size L of the diffusive halo of the Galaxy, 
ii) the normalization of the diffusion coefficient Ko, iii) its 
slope 5 and iv) the constant galactic wind Vc. The boundary 
conditions are N{z = L,r) = and N{r ^ R, z) ~ 0, where 



Set 


S 


A'o {kpc2 Myr-i) 


L (kpc) 


Vc (km s-l) 


max 


0.46 


0.0765 


15 


5 


best 


0.7 


0.0112 


4 


12 


min 


0.85 


0.0016 


1 


13.5 



Tabic 1 

Prop agation parameters consistent with B/C data dMaurin et al.L 
I2OOII) . The set labeled best corresponds to the best fit to B/C data, 
while those labeled min and max co rrespond to sets wh ich give 
minimum and maximum exotic fluxes l lDonato et al.L |2004|) . 



tion, as given, e.g., in Berezinskii et al. ( 199GI ). Through- S 
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Fig. 1. Reference interstellar standard p flux along with its fit function 
as given by Eq. (2). 

the Galactic radius is i? = 20 kpc. The propagation param- 
eters h ave been deter mined from the analysis of the B/C 
ratio in lMaurin et al.l (2001|): only a sub-region of the prop- 
agation parameter space is compatible with B/C data. The 
accuracy of the approximate formulae given in this work is 
studied for this sub-region only. Within the latter, the sets 
of parameters lea ding to the maximum /minimum primary 
exotic fluxes (see Donato et al.l 12004 and the companion 
paper) are gathered in Tab. 1, along with the configuration 
corresponding to the be st fit to B/C data. These three sets 
were previously used in iDonato et al. ( 2004 ) ; Barrau et al 



(200. 



2.3. Secondary p fl,ux 

The secondary p flux has been studied in detail in 
Donato et al.l (j200l[ ) : every set of propagation parameters 
consistent with the observed B/C ratio leads to the same p 
secondary flux, within the ^ 10% uncertainties mostly of 
nuclear origin. Moreover, the p fluxes o btained from other 
mode ls are in fair agreement (see e.g. iMoskalenko et al 
2OO2I) . 
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If one wishes to exclude, for example, an exotic primary 
contribution by scanning the SUSY parameter space, it 
is useless to calculate again and again the same quantity 
with various input propagation parameters. The following 
fit function, once modulated, provides a proper description 
of the measured p secondary flux, regardless of the propa- 



gation model. Setting x = \og{Ek) where Ek is the kinetic 
energy in GeV, the standard IS p flux in (GeV s sr)"-'^ 
is parameterized as 



standard 



exp C.x' j if Ef, < 11 GeV 
exp [Dqx^^ ) otherwise ; 



(2) 



The constants are defined by Co = -3.211, Ci = 0.12145, 
C2 = -0.2728, C3 = -0.075265, C4 = -0.007162, and 
-Do = -2.02735 and Dl = 1.16463. The fit is shown in 
Fig. 1. 

2.4. Exotic primary p flux 

The derivation for the primary flux has been presented 



Barrau et ah (|2002P1 . The solution can be found using 



R 



a Bessel expansion (see e.g.. lMaurin et al.ll200l[ ) 

7VP■^'-(r,z) = f;A/^P"-(z)Jo(c 
where Ci is the i-th zero of Jo and 

A/;P-(.)^M(0)e^ -^""jyL7)] 

smh (^) 

v^(L^,) Smh{S^z/2) y,{L) _ yjjz) 
sinh(S',L/2) KS, KS, 



(3) 



Each Mi (0) obeys a differential equation 
2h 



M(0)=A/^*(0) 



(4) 



dE ' ' dE^ 
where A/"* (0) is the solution without energy redistributions 

-VaL/2K yi{L) 



A,sinh{S^L/2) 



(5) 



E quation (4) is s olved with a numerical scheme detailed 



in iDonato et al.l (|200lr ). The tertiary contribution is not 
shown in the former equation although it is taken into ac- 
count in our numerical calculations. The remaining quan- 
tities are defined as follows: 



2 / e 
/o 



2K ^'^ ^ 



sinh 



Si . I , 



Ai=2hTtot + Vc + KS^ coth 



qi{z')dz' 



In the above expression, qi{z) are the Fourier-Bessel coeffi- 
cients of the source term iSprim(^, z). Some numerical issues 
about the convergence of the series, when using cuspy dark 
matter halo, are underlined in App. B. 



Fo otnote that there is a misprint in Eq. (A. 5) of l l Barrau et al.l . 
I2OO2I) . The product KAiSi should read KSi. 



Halo model 


a 


/3 


7 


Tc [kpc] 


Cored isothermal 


2 


2 





4 


NFW CNavarro et al.. 1997. 2004) 


1 


3 


1 


25 


Moore (Moore et al., Diemand et al.) 


1 


3 


1.2 


30 



Table 2 

Parameters to be plugged into Eq. (7) to obtain various modelings 
of the dark matter distribution profile in the Milky Way. 

2.5. Dark matter profile 

The primary exotic source term follows the generic form 

5prim('', 2) = Q^iE) X /source z) , (6) 

where the spatial dependence of the source term is normal- 
ized to 1 at Solar position. For evaporation of Primordial 
Black Holes or for SUSY-like particles, /sourcc(?', z) is equal 
to 

/Dark(?',2) or /Dark(''>2;) , 

where fuiukir, z) is the dark matter distribution profile, 
also normalized to 1 at Solar position. Any information such 
as, e.g., the relic density, the fraction of the dark matter 
halo filled by the new candidate or any quantity related to 
the dark matter candidate is irrelevant at this stage: it is 
absorbed in the energy dependent term Q{E). The right 
factors should be implemented accordingly when using the 
approximate formulae. 

The following generic form is taken for the profile: 

, (/3-7)/a 

Rq 



fDm-kir,z) 



^/r"^ -\- z"^ 



Rq 



rj? + \/r^ + z^ 



(7) 

where parameters are given in Tab. 2. The distance to the 
galactic center is set to Rq = 7.5 kpc, a value that s everal 
methods seem to converge to (jNishivama et al.l . l2006l ). This 
is at variance with th e usually recommended Rq = 8.0 kpc 
(e.g. Yao et al. 20061 ). but has no impact on the derived 
results below. Indeed, Rq is just another parameter, and 
our approximate formulae will perform well for any user- 
preferred value « 8.0 kpc. 

3. Approximate formulae 

The previous formulae (see Sec. 2.4) can be tricky to im- 
plement in practice (numerical inversion, tertiary contri- 
bution, convergence issues), especially when energy redis- 
tributions are included. In this section, energy gains and 
losses are discarded and a simplified formalism is presented 
for the calculation of exotic primary p fluxes. Some of the 
formulae given hereafter are well suited for extensive com- 
putation. 

3.1. Propagator (no side boundaries: R oo) 



In lTaillet et al.l lj2004l ) . the expression for the propagator 
corresponding to the same 2D diffusion equation (without 



Q 



energy redistributions) was extraeted in the limitT^ R oo. 
Henee, the problem is invariant under radial translations 
and we choose the location of the observer as the origin 
(r = 0). The propagator Q{r ~ r') is solution of 

d 1 

-A'A + + 2hnS{z)v(7 > G ^ S{r - r') 

and reads (see App. A) 



g{r, z) 



(8) 



2ttKL 

oo 

c^^KQ{ry^k^^ + fc,2) sinfc„i sinfc„(L — z) 

where Kq is the modified Bessel function of the second kind, 
fc„ is solution of 



2fc„ cos kjiL = — fcfj sin fc„L 



Cn is defined as 



and 



sin knL cos fc„i 
knL 



(9) 
(10) 



fc, = Vc/{2K) . 

The quantity k^ depends on the total destruction rate Ftot 
and the constant convective wind Vr-: 



kd = 2h Ttot/K + 2k,, . 
The flux at solar location is given by 



(11) 



Nq = 2 / / g{r,e,z)Sp,in,{r,e,z)rdedrdz . 

Jq Jo Jo 

where 5prim('"j 9, z) is the dark matter source term as seen 
from the Solar System. Hence, a single (numerical) integra- 
tion is needed in the propagator approach: this represents 
a great benefit compared to the full 2D description given 
in Sec. 2. 



3.2. The ID model 

Another way to tackle the simplified problem (neglect- 
ing energy redistributions) is the use of a ID model. We 
now consider an infinite disk whose density and source dis- 
tribution do not depend on r. The constant wind diffusion 
transport equation for a generic source term Q(z, E) = 
q{z)Q{E) without energy redistributions reads 

d'^ N dN 

2h5{z)TtotN ^ q{z)Q{E) . (12) 



K- 



dz^ ' ^ dz 

As we do not take into account energy losses, the energy 
E is decoupled from the spatial dependence z so that it 



^ The propagator f or fi nite R was presented and used 
in |m aurin and TaiUetl ||2003|) to obtain the most Ukely origin of ex- 
otic p detected on Earth. This propagator relied on Bessel expansion. 
For practical use, i? 3> 1 was required, so that it was equivalent to 
the form presented below. However, in the context of a quick and 
simple formulae, it is largely outperformed (in term of convergence 
properties) by the one presented in this paper. 
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Fig. 2. From top to bottom, the max, best and the min propagation 
sets (sec Tab. 1). The dashed lines correspond to the propagated 
fluxes without energy redistributions, while solid lines correspond to 
the same propagated exotic fluxes, but with energy redistributions 
(energy losses, reacceloration, tertiary source term) taken into ac- 
count. The shaded regions correspond to p which cannot reach Earth 
(because of solar modulation, see text). The two panels correspond 
to different dark matter profiles and noutralino masses. 

is omitted. In ID-models, only z-dependence is allowed. 
For simplicity, it is further assumed that the source term 
q is constant throughout the Galaxy. We refer the reader 
to the companion paper for details of the calculation. The 
analytical solution reads 



qL \{l + a + £,)■{! 



Vc I a + e(l-e-") 

qL r 1 - (1 + a)e- 



z 
L 



(13) 



Niz) 

where ^ = hVtotL/K and a = VcL/K. 
4. Results 



In this section, the validity of the previous approximate 
formulae is checked. We first study the influence of en- 
ergy redistributions within the 2D numerical framework 
and then compare the two simplified formulae to the 2D 
numerical model (without energy redistributions). 

4.1. Influence of energy redistributions in the 2D model 

The realistic dark matter profile (as described in Sec. 2.5) 
is implemented in the full 2D-model (see Sec. 2.4), with 
and without energy redistributions. When energy redistri- 
butions are neglected, the quantity A/i(0) to plug in Eq. (3) 
becomes the quantity A/7(0) defined in Eq. (5). We use the 
generic term error for the difference between solutions cal- 
culated with and without energy redistributions. 



A 



4.1.1. Expected effect of energy redistributions 

As energy redistributions act on the spectrum, it is 
not obvious to disentangle which of the following in- 
put parameters — source spectrum, profile, transport 
coefficients — is the most crucial regarding error. On a gen- 
eral footing, energy losses and tertiary contributions tend 
to replenish the low energy tail in case of a dropping flux at 
these same energies (this has a partic ularly important im - 
pact on the standard secondary p, see lDonato et al.ll200l[ ). 
On the other hand, reacceleration tends to smear the spec- 
trum at low and high energy in case of a flux peaking at 
GcV energy. Taking all these points into account, the fol- 
lowing remarks — regarding energy redistributions — can be 
made: 

Source spectrum: in most cases, PBH-like source 
sp ectra do not decre ase at low energy (e.g.. Fig. 3 
of Barrau et al. I I2OO2I) . whereas SUSY-like candidates 



Isothermal profile 



Moore profile 



are quite flat (e.g.. Fig. 6 of lDonato et al.ll2004l ). Hence 



taking various source spectra is not expected to play a 
major role in error. 

Dark matter profile: changing the spatial distribution 
does not affect the spectral distribution: it only changes 
the normalization of the propagated spectrum. Hence, 
no significant effect on error is expected when switching 
between SUSY-like to PBH-like source terms. 

Propagation parameters: Large winds dominate over 
diffusion as the energy diminishes: they decrease low en- 
ergy fluxes, making them more sensitive to energy redis- 
tribution effects. These large wind happen to be asso- 
ciated with a small halo size and a small reacceleration 
parameter. This effect is larger than the two others (see 
below) . 



4.1.2. Behavior of error 

Figure 2 displays the two resulting fluxes for various in- 
put configurations. The decrease of the flux observed at low 
energy for the best and minimal propagation sets is due 
to the Galactic wind. For the smallest L, the propagated 
flux is steeper (because of larger wind values, see Tab. 1), 
hence more affected by energy redistributions. Notice that 
only for L ~ \b kpc (corresponding to = 5 km s~^), can 
the shape of the source spectrum be truly appreciated. As 
expected, the error calculated when varying the ncutralino 
mass or the dark matter profile is less important. 

The error can be as large as a factor of '--^ 5 at very low 
energy, but we remind that these are p at IS energies, which 
cannot reach Earth because of solar modulation (shaded 
region in Fig. 2). Indeed, the top of atmosphere energy 
Etoa is related to the interstellar energy Eis by Eiq = 
Etoa + wliere^ is th e modulation parameter e [500 — 
1000] MV (|Perki[l983). 

Above a few hundreds of McV (lowest IS energy reach- 
able for low solar activity), the difference obtained is al- 
ways < 50%. Above a few tens of GeV, energy redistribu- 
tions are always negligible, so that error tends to zero. In 
between, the only situation when error can be large corre- 
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Fig. 3. Both figures correspond to the ratio of the p exotic flux to its 
source term (in arbitrary units) as a function of the kinetic energy. 
Lines: 2D-semi-analytical model, no energy redistributions (Sec. 2.4) 
for PBH-like (solid lines) or SUSY-like (dashed lines) sources. Sym- 
bols: Eq. (13), ID-model (independent of the dark matter profile). 
The three sets of propagation parameters are those of Tab. 1. Left 
panel: Isothermal profile; right panel: Moore profile. 

sponds to a steep spectrum dropping rapidly to zero, which 
is uninteresting in the perspective of detecting any excess. 

4.1.3. Summary 

The model neglecting energy redistributions is accurate 
enough to determine exotic fluxes. This result is in agree- 
ment with the intuitive idea that exotic species, in order 
to reach us, cross the gaseous disk less often that standard 
ones, implying that they are less subject to energy redistri- 
butions. 

4.2. Approximate formulae vs 2D reference model 

4.2.1. Propagator 

As expected, the propagator (see Sec. 3.1) and the 
Bessel expansion (see Sec. 2.4) formalisms are consistent, 
except when the halo size has the same extension of the 
radial extension of the Galaxy R. In that case, we no 
longer have L R and the si de boundary starts affect- 



ing 



the flux. 



Maurin and Taillet 



as di scusse d 
2003). 



in (Taillet and Maurin. 2003 



4.2.2. ID-model 

Figure 3 compares the ID-model (at z = 0) to the 2D 
semi-analytical model: 

- For PBH-like sources, as long as L < 7 kpc, the ID- 
model provides a good description of exotic p fluxes. 

- For SUSY-like sources, ID-models are valid for very 
small halo size L only. 

This good agreement is easily understood if we remember 
that the closest boundary defines a cut-off distance: anti- 
nuclei from sources located further away are exponentially 



supp ressed (ITaillet and Maurin . l2003t iMaurin and Tailletl . 
20031 ). Hence, we are only sensitive to the average value of 
/Dark withiu this cut-off distance, close to one at solar po- 
sition, as taken in the ID-model. As soon as contributions 
close to the Galactic center are less suppressed, this equiva- 
lence breaks down; all the more for SUSY-likc source terms 
compared to PBH-like source terms (because of the square 

in /Dark)- 

5. Numerical routines 

Taking advantage of the two simplified formalisms de- 
picted above (propagator and ID model), routines are pro- 
videc{3, to evaluate the p flux from any input dark matter 
profile, source spectrum, and propagation parameters. The 
main routine returns the ratio of the propagated IS exotic 
flux (for a user-specified source term) to the standard p IS 
flux, i.e. 

exotic/ standard ' 

The input parameters are the following: 

- Propagation parameters: (kpc^ Myr~^), S, Vc 
(km s~^) and L (kpc); 

- Exotic source spectrum: Qoxoticl^fc) (GeV~^ m~^ s^^) 
and the corresponding kinetic energy Ek (GeV); 

- Dark matter profile: the parameters a, /3 and 7 (see 2) 
as well as the source type (SUSY-like or PBH-like). 

The resulting fluxes, in (GeV m^ s sr)^^, are obtained at 
solar position (i?Q , 0) . The IS background p spectrum is 
described using a fit function (see Sec. 2.3). We chose to ex- 
press all fluxes as IS quantities to free the user from using 
a modulation routine. There are uncertainties associated 
with the choice of the modulation/demodulation scheme, 
but these are not la rger than the error mad e using approxi- 
mate formulae (see iMoskalenko et al"]l2002l for a discussion 
on solar modulation). 

We remind that the range of propagation parameters 
given in Tab. 1 implies that large uncertainties remain on 
exotic fluxes (without mentioning the fact that the constant 
wind model is probably not the definitive model for the 
Galaxy). To make conservative estimates, we recommend 
the reader to use the min set of propagation parameters 
associated with the condition O > 1 to exclude an exotic 
model. Regarding the most likely p exotic flux in this model, 
it is given by the best set (L — 4 kpc, sec Tab. 1). In any 
case, it has to be kept in mind that energy smaller than 
~ 600 MeV IS correspond to p which arc almost never 
detected on Earth. 

6. Conclusion 

Given the present accuracy of propagation parameters, 
exotic anti-proton fluxes suffer large uncertainties (see also 
the companion paper). This limits the benefit of a detailed 



calculation, which can be expensive in term of computa- 
tional power, when repeated thousands of times. 

We showed that discarding energy redistributions 
(coulomb, ionization, adiabatic losses, reacceleration and 
tertiary contribution) provides exotic p fluxes accurate at 
the level of < 50% for IS energies greater than ~ 1 GeV 
(leading to > 100 MeV energies once modulated). This 
opens the possibility to use one of the following simpler 
formulae: 

- the approximate ID-model, which gives the correct flux 
as long 'AS L <7 kpc for PBH-like sources, or L < 1 kpc 
for SUSY-like sources; 

- the 2D-propagator. 

The advantage of the ID-modcl is its compactness and sim- 
plicity (no dependence on the dark matter proflle), while 
the second approach is more accurate. The former could be 
preferred in situations where speed is a decisive factor, e.g. 
when scanning a large parameter space as in SUSY studies. 
Moreover, in order to be conservative in exclusion studies, 
small halo size should be taken. In particular conservative 
estimates are made when using the 3''*^ line of Tab. 1: this is 
exactly the regime where ID-formulae hold. We also sup- 
ply numerical routines that could be easily implemented in 
indirect detection codes. 

So far, these formulae arc valid only in the framework of 
the constant wind model. Although all the discussions were 
developed with respect to exotic p fluxes, all the results 
can be transposed almost without any modiflcation to the 
d case (the total cross section has to be taken accordingly 
to the species considered). 
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Appendix A. Check of propagator formulae 

The validity of the propagator is checked in a simple case, 
by integrating it over a constant source term q(r, z) = 1 
(ID formula). 

A.l. Purely diffusive regime: 14 = 

A. 1.1. No spallations, kd ~ 

In that case, fc„ = (n + l/2)7r. Integration over z (using 
g(r, z) = 1) leads to 

dzg{r,z)^-—Y,^-j^Ko{knr), (A.l) 

n=0 " 

and integration over the disk 2-Kr dr gives 



n=0 " 
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yKa{y)dy 



a 



Using the properties J y Ko{y) dy 
yKi{y) 1 when y 0, 

i-ir 



-yKi{y) and 



ri=0 " 



Finally, for i? ^ 1, the expression becomes 



n=0 



(-1)' 



2L 



2 °o 



(-1)" 



7r3A- (n + 1/2)3 



The infinite sum has the value + 1/2)^ = 

7r3/4 « 7.7515, so that 



ID — model 



(0). 



Note that, in the infinite sum, the first term n = is 
equal to 8, so that it only slightly overshoots the exact 
value. Hence, only a few terms are required in the sum to 
converge quickly to the correct value. 

A. 1.2. With spallations, kd = 2hTtot/K 

There is no longer a simple expression for fc„. Following 
a similar derivation as the one above, we find 



iV(0) 



4 / sin fc„ L cos fc„ L 

El-- 



KkdL 



kn.L 



cos kri L 



[cos knL — 1] 



(A.2) 



A numerical check of this sum confirms that Eq. (A.2) 
equals the ID-model, i.e. Eq. (13) where Vc = 0. This 
means that 

cos{knL) [cos(fc„L) — 1] ktjL 



^ knL {knL - sin(A:„L) cos(fc„L)) 4(2 + kdL) 



A.2. Diffusive/ convective regime: Vc /= 

The integration leads to quite similar results as for the 
previous case (with spallations). A numerical check con- 
firms that Eq. (13) is recovered. 

Note that the integration should be performed for sources 
located at all r, with r — > oo. In practice, it is sufficient 
(and it saves a lot of computational time) to integrate only 
from to lOL (or from to min(10i, lOicff)) in the case of 
a convective wind, where L^ff = K/Vc- In that case, a grid 
of ~ 25 points for all coordinates (r, z and 9) is sufficient 
to reach the correct result. 

Appendix B. Numerical instabilities in the 2D 
model 

Dealing with Bessel functions is a source of numerical in- 
stabilities. These occur if i) i? ^ L or ii) too many Bessel 
orders are used, or even, in the case of cuspy dark matter 
halo, if too few orders arc used. As the propagator formu- 
lation allows cross-checks, we take the opportunity to give 
a few recommendations regarding the parameters to use: 



(i) always keep R in the range 5L* < i? < lOL*, with 
L* = mm{L,K/Vc). The limit lOL"^ is set to avoid 
numerical instabilities w hile the limit 5L* aUows to 
avoid boundary effects ( Maurin and Taillet . 20031 ) 
(the boundary in R decreases the fiux at most to a 
few tens of percent for a Moore profile). 

(ii) For a smooth profile (e.g. Isothermal), 20 Bessel func- 
tions give an excellent convergence. 

(ui) For NFW and Moore profiles in the case of SUSY- 
like candidate (50 Bessel functions are sufficient for 
P BH-like candid ate), following a procedure inspired 
bv lBarrau et al.l f|2005, ) . we replace in the calculation 
the standard profile, by 



5Dark('') = 



'th 



/Dark('"th) • TT T 



2 sin(7rr/rth) 



nr/rth 



if r < rth ; 

(B.l) 



f2 

J Dark 



otherwise; 



with T = (3 — 27) ^ and 7 is the slope from Tab. 2. In 
these expressions, r = 



''cyi + -^^ denotes the spheri- 
cal coordinate. In any case, we find that, setting rth = 
400 pc, 50 Bessel functions give a good accuracy, while 
100 functions allow to reach an excellent convergence. 
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